Draft version February 12, 2013 

Preprint typeset using I^'T^]X style emulatoapj v. 5/2/11 



o 

(N 



O 



(N 

a^ 

(N 

(N 
O 
m 



AN INSTABILITY DUE TO THE NONLINEAR COUPLING OF p-MODES TO 5-MODES: IMPLICATIONS FOR 

COALESCING NEUTRON STAR BINARIES 

Nevin N. Weinberg^ Phil Arras^, and Joshua Burkart^ 

^Department of Physics, and Kavli Institute for Astropliysics and Space Research, Massachusetts Institute of Technology, Cambridge, 

MA 02139, USA; nevinOmit.edu 
^Department of Astronomy, University of Virginia, P.O. Box 400325, Charlottesville, VA 22904-4325, USA and 
^Department of Physics, 366 LeConte Hall, University of California, Berkeley, CA 94720, USA 

Draft version February 12, 2013 

ABSTRACT 

A weakly nonlinear fluid wave propagating within a star can be unstable to three-wave interactions. 
The resonant parametric instability is a well-known form of three- wave interaction in which a primary 
wave of frequency uJa excites a pair of secondary waves of frequency ujb + lOc — oJa- Here we consider 
a nonresonant form of three-wave interaction in which a low-frequency primary wave excites a high- 
frequency p-mode and a low- frequency f;-mode such that wt, -|- Wc 3> Wa • We show that a p-mode can 
couple so strongly to a g-mode of similar radial wavelength that this type of nonresonant interaction 
is unstable even if the primary wave amplitude is small. As an application, we analyze the stability of 
the tide in coalescing neutron star binaries to p-g mode coupling. We find that the equilibrium tide 
and dynamical tide are both p-g unstable at gravitational wave frequencies /gw ^ 20 Hz and drive 
short wavelength p-g mode pairs to significant energies on very short timescales (much less than the 
orbital decay time due to gravitational radiation). Resonant parametric coupling to the tide is, by 
contrast, either stable or drives modes at a much smaller rate. We do not solve for the saturation of 
the p-g instability and therefore we cannot say precisely how it infiuences the evolution of neutron 
star binaries. However, we show that if even a single daughter mode saturates near its wave breaking 
amplitude, the p-g instability of the equilibrium tide will: (i) induce significant orbital phase errors 
(A(/> > 1 radian) that accumulate primarily at low frequencies (/gw ^ 50 Hz) and (ii) heat the neutron 
star core to a temperature of T ~ 10^*^ K. Since there are at least ^ 100 unstable p-g daughter pairs, 
A(j) and T are potentially much larger than these values. Tides might therefore significantly infiuence 
the gravitational wave signal and electromagnetic emission from coalescing neutron star binaries at 
much larger orbital separations than previously thought. 

Subject headings: binaries: close - hydrodynamics - gravitation - stars: interiors - stars: neutron - 
stars: oscillations - waves 



1. INTRODUCTION 

In the standard theory of stellar oscillations, one 
solves the linearized fiuid equations with the assump- 
tion that the waves that propagate within a star do not 
interact with each other. However, in some stars in- 
ternal waves are driven to such large amplitudes that 
the linear approximation becomes invalid. Examples 
of systems in which nonlinear wave interactions can be 
important include close binaries (iK uniar & Goodman 



1996; Goodman & Dickson"1998';13arkcr & Oeilvic 201C 
Weinberg ct al. 2012; Fuller & Lai 2012; Burkart et aj] 
2012} ). the su n (Kumar fc GoldreichI 1 1 9891). white 
dwar fs (|Wu fc Goldreich..200"l . RR Lvrae (|Molnar et al.l 
I20121 ). and neutron stars that are either newly 
formed (iW einberg fc Ouataertll2008h or rapidly rotating 
()Arras eta l. 2003). 

As long as the nonlinearities are not too strong, the 
wave interactions can be described using a perturbative 
approach. At the lowest nonlinear order, the interac- 
tions involve three-wave couplings in which a large am- 
plitude parent wave a excites pairs of daughter waves 
(5, c) . These interactions often occur as a resonant para- 
metric instability in which the parent's (possibly driven) 
oscillation frequency coa nearly equals the sum of the 
daughters' iiatura l oscillation frequency uib-\-uJc (see e.g., 
iHasselmannl 119671 ). The studies cited above exclusively 



considered such resonant interactions between parents 
and daughters. Here we consider the stability of non- 
resonant interactions among strongly coupled waves. 

The nonlinear coupling strength is sensitive to the spa- 
tial structure of the waves. For three-wave interactions, 
the coupling strength is parametrized by the coupling co- 
efficient Kabc- The magnitude of Kabc is largest in regions 
where the radial wavenumbers satisfy |fcf, — /cd ^ ^a, the 
usual condition for momentum conservation; otherwise, 
the waves are incoherent and t heir spatial oscillations 
tend to cancel out the interaction (|Wu fc Goldreichl[2001l : 
IWeinberg etHllml hereafter WAQB). Thus, 

Kabc can 

be large for even long wavelength parents if kt ^ kc- For 
waves of high radial order and low angular degree, the 
dispersion relation is fc ~ w /cg for high-frequency acous- 
tic waves (p-modes) and fc ~ AN/ur for low-frequency 
internal gravity waves (g-modes; here Cs is the adiabatic 
sound speed, A^ = £{£ -\- 1), i is the angular degree, TV is 
the Brunt- Vaisala frequency, and r is the radial coordi- 
nate). There are therefore three ways to satisfy kb ~ fcc 
at a given radius: (i) the daughters are both g-modes 
and AcLOb — AbWc, (ii) the daughters are both p- modes 
and LOb — i^c, or (iii) one daughter is a p-mode and the 
other is a g-mode and 
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where we took b to be the p-mode and c to be the g- 
mode. In cases (i) and (ii), the three waves can satisfy 
the nonhnear resonance condition if ui, — ujc — Wa/2. In 
case (iii), ujb ^ Wc and the resonance condition cannot 
be satisfied (assuming that < uJb', e.g., the parent is a 
(7-mode or hnearly driven by a tide). We will show that 
p-modes can couple so strongly to .g-modes that such 
nonresonant interactions can, nonetheless, be unstable 
even for relatively small amplitude parent waves. 

The primary application of p-g mode coupling that 
we consider in this paper is in the context of tides in 
coalescing neutron star-neutron star (NS-NS) and neu- 
tron star-black hole (NS-BH) binaries. These binaries 
are the most promising sources for ground-based gravi- 
tational wave o bservatories such as LIGO and VIRGO 
(jCutler fc Tho rnc 2002). Within the next few years, 
advanced versions of these detectors should be taking 
data with sufficient sensitivity that they will detect the 
first gravita tional wave signatu re of a compact-binary 
coalescence ()Abadie et al.|[20To() . Extracting the signal 
from detector noise requires accurate theoretical tem- 
plates and thus a precise understanding of the gravita- 
tional waveform. If a phase error of even « 1 radian ac- 
cumulates over the final w 10^ or bits, it can significa ntly 
decrease the source detectability ([Cutler et al.lll993l) . 

Tidal interactions extract energy (and angular momen- 
tum) from the orbit and, depending on the nature and 
rate of the internal dissipation, deposit it within the star 
as some combination of mode and thermal energy (and 
mode and spin angular momentum). As a result, tidal 
interactions in NS binaries modify the rate of inspiral 
and lead to phase shifts in the gravitational waveform 
that may affect sour ce detectability if sufficiently large 
and unaccounted for (jBildsten fc Cutleilll992t iKochaneM 
IT992) . Conversely, tide-induced phase shifts may encode 
highly sought information about the NS equation o f state 
(jPlanagan fc Hindereiir2008HHinderer et al.ll2010[ ). 

As the NS inspirals, tides induce a large scale dis- 
tortion of the star (referred to as the equilibrium tide) 
and also excite resonant oscillation modes (the dynam- 
ical tide). Although the equilibrium tide stores consid- 
erable energy, if its dissipation is determined entirely by 
line ar processes, it does not sign i ficant ly affect the or- 
bit (|Bildsten fc Cutleii [1991 iLail [1991 . The dynami- 
cal tide was first studied in non-rotating NSs, where the 
resonant mod es are q-m odcs with frequency < 100 Hz 
((Rciscncggcr fc Goldrei ch 1994; Lai 1994). These stud- 
ies found that while the effect on the gravitational wave- 
form is small, linear dissipation of the excited g-modes 
can heat the NS core to a temperature of ~ 10* K. Fur- 
thermore, rapid rotation can strongly enhance the tidal 
effects and lead to the excitation of r-modes and inertial 



wav es, resulting m p nasc smtts 01 ~ U.i to 3> i radi- 
ans (iHo fc Lail[l999t[Lai fc Wu 2001 [Flanagan fc Racind 



12007 ). However, the required spin frequencies are higher 
than is thought to be likely for NS-NS binaries. 

All of these studies ignored nonlinear interactions and 
instead assumed that linear theory is valid at gravita- 
tional wave frequencies /gw < 400 Hz [3 Some of the 



^ While several groups now carry out hydrodynamic numerical 
simulations of compact object inspiral using realistic e quations of 
state (e.g.. [DecRslin et aLlfeooTI : ISekiguchi et al.ll2011l) . they only 
simulate the last few orbits before the merger (/gw 400 Hz). 



studies argue that because the amplitude of the tidal per- 
turbations are ^ 1% of the NS radius at these frequen- 
cies, the linear approximation should be valid. However, 
the validity of the linear approximation depends on more 
than just the amplitude of the perturbations; it also de- 
pends on the strength of the nonlinear coupling Kate be- 
tween the primary perturbation and other modes within 
the star, and the individual properties of those modes 
such as their frequency and linear damping rates. As a 
result, even if the perturbation amplitude is <i; 1% of 
the stellar radius, it is potentially unstable to nonlinear 
wave interactions. Indeed, we will show that by the time 
a binary first enters LIGO's bandpass (/gw ~ 20 Hz), 
the equilibrium and dynamical tides are unstable to non- 
linear three-wave interactions. These instabilities drive 
rapidly growing, short wavelength modes and can po- 
tentially lead to significantly enhanced tidal dissipation 
relative to linear theory predictions. 

The structure of the paper is as follows. In §[2]we carry 
out a stability analysis of p-g mode coupling for a linearly 
driven parent wave. In § [3| we calculate the strength of 
p-g mode coupling in NSs. In § [3| we show that coa- 
lescing NS binaries are subject to the p-g mode coupling 
instability (PGI), with the equilibrium tide and dynami- 
cal tide serving as parent waves that drive the daughters. 
We also evaluate the stability of the tide to the resonant 
parametric instability. In § [5| we estimate how the PGI 
might influence the orbit and tidal heating of coalescing 
NS binaries. We conclude in § [B| and briefly discuss the 
potential influence of the PGI in systems other than NS 
binaries. 

2. STABILITY OF p-g MODE COUPLING 

Consider a linearly driven parent a that is coupled to a 
daughter pair (5, c) . To lowest nonlinear o rder, the am- 
plitud e equations take the form (see, e.g.. ISchenk et al.l 
[200I and WAQB) 



9a - 
96 



■ laQa + OJala^i^a [Ua{t) + l^*abcqU*c] (2) 
- -fbClb + t^lqb = ^l^lbSlll (3) 



9c + 7c9c + i^lqc = ojlKl^^q^Qb , (4) 

where the q's are complex mode amplitudes, the 7's are 
linear damping rates, and we assume that the parent is 
harmonically driven as Ua{t) — J/aC""^'. In order to 
determine whether the parent is stable at its linear am- 
phtude, let qb,c = Ab^cC^""'"*'^''''''" where s is complex and 
<7b.c is real. If ab <Jc = ^, the time dependences cancel 
in the amplitude equations and we obtain the character- 
istic equation for s. For daughters that are not resonant 
with the parent {ujb ^ Wcw), the solution to the char- 
acteristic equation states that if (see Appendix [^ 



\qa\ > \Kabcl 



(5) 



the daughters are unstable and grow exponentially at a 
rate 

Tbc^UJclKabcqal- (6) 

The stability criterion applies even if the linear damp- 
ing rates are large, i.e., jb,c ^ ^b.c- Furthermore, even a 
static (cj = 0) parent can be unstable(3 Such an insta- 

^ The instability of an oj = parent implies that the static 
background is unstable to small perturbations (as in a Rayleigh- 
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bility has the features of what iWul ()1998l ) refer to as an 
ampUtude instabihty. 

Here the ampUtude equations are second order in time 
whereas in WAQB they are first order in time. This is 
because here we adopt a configuration space mode expan- 
sion of the form ^ = Qa^ai where ^ is the Lagrangian 
fluid displacement and the sum is over all modes. WAQB 
instead adopt a phase space mode expansion of the form 

= J2a1a.i^a^~'^^a^a) where the sum is over aU 
modes and their complex conjugate. Both forms of the 
amplitude equations are, of cours e, valid and describe the 
same physics (see Appendix C in lSchenk et al.l[2002t ). 

By contrast with the non-resonant case, resonant 
daughters (wb + cUc — oj) are unstable if 



1 9a 



> 



Iblc 



Kobe Wfc'^c 



A2 



-I 1/2 



(76 + 7c) 



(7) 



where A 



Wb -t- UJc 



u! is the detuning (see, e.g. 



INishikawal[l968[ WAQB). If |A| < 7b,c < i^b^c, the fac 
tor multiplying |ko6c|~^ in the above equation is <C 1. 
Nonetheless, if the non-resonant Kate is much larger than 
the resonant Kate, the non- resonant instability can have a 
lower amplitude threshold. This is indeed the case in coa- 
lescing NS binaries: in § I3.2l we show that the equilibrium 
tide couples to a (p, g) daughter pair with 10^ times 
greater strength than it couples to a resonant daughter 
pair (independent of orbital frequency). As a result, the 
equilibrium tide is unstable to the PGI but not the res- 
onant parametric instability. 

Note that the resonant instability criterion derived in 
WAQB (eq. [7]) does not reduce to the non- resonant 
criterion (eq. [5]) in the limit A — > Wb. In deriving the 
resonant criterion, WAQB only accounted for the sum 
over modes and not their complex conjugate; as a result, 
equation ([7]) is only accurate for nearly resonant daughter 
modes Wb + ~ w. 

We illustrate the PGI in Figure [T] by numerically solv- 
ing equations ©-([l]) for a non-resonant, unstable three- 
wave system. The modes are each given different ini- 
tial conditions (amplitudes, phases, etc.). We find that 
the instability is not sensitive to the initial conditions 
and that the daughters quickly attain similar amplitudes 
\qh\ ~ \qc\ once the nonlinear coupling terms dominate 
the linear terms. 

Near the instability threshold \qa\ ^ \nabc\~^ and 
the linear terms are comparable in magnitude to the 
second-order nonlinear terms (left and right hand side 
of eqs. respectively). One might therefore won- 

der whether higher-order terms, which we neglect, are in 
fact important near threshold (or even whether pertur- 
bation theory is breaking down). We show in § 13.41 that 
this is not the case. 

3. STRENGTH OF p-g MODE COUPLING 

The possible types of p-g couplings include those with 
(p, g) daughter pairs 

9-P9, 



p:pg, tide:pg, 



Taylor instability). Note, however, that we do not consider lu = 
parents in this paper; even the equilibrium tide is time-dependent 
(e.g., the £ = 2,m = ±2 harmonic oscillates at /gw). 
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Fig. 1. — Amplitude evolution of a three- wave system in which a 
parent mode a (upper solid line) is linearly driven at a frequency 
ui and is coupled to a pair of daughters (b, c), where 6 is a p-mode 
(dashed line) and c is a g-mode (lower solid line). Time is shown 
in units of 2it/u}. The parameters of the system are K^bc = 10^, 
Ua = 3 X 10"^, and in units of u), {(^a,^a) = (1.005,0.01), 
(ujb,-rb) = (104,1.0), (iJc,7c) = (0.7,0.3). The modes are given 
initial amplitudes of log \q\ = —10, —14, —12, respectively. As the 
parent approaches its linear amplitude of 2 X lO"'^, it crosses the 
amplitude (dotted line), and the daughters, which initially 

decay, become unstable and grow at a rate Fbc — i^clnabcQal- 

and those with {p,p) or {g,g) daughter pairs 
g:pp, p:gg, 

where a : be denotes the coupling of a parent a to a 
daughter pair (6, c) with modes types p, or a tide 
(equilibrium or dynamical; we also show results for /- 
mode parents). First consider the couplings involving 
{p, g) daughter pairs. As described in § [1] the constraint 
1^6 ~ I ^5 implies that their frequencies must satisfy 
UbUJc — ^cNcs/r. The spatial extent of the coupling re- 
gion, and therefore the magnitude of Kabc^ depends on 
how rapidly Ncg/r varies with radius. In stellar cores, 
often N (X r and Cs ~ constant and the modes couple 
well throughout the core. The coupling can be partic- 
ularly strong in stars with large cores, such as NSs and 
white dwarfs. In Figure [2] we show the structure of an 
M = 1.4Mq (i? ~ 11.7 km) NS assuming the Skyrme 
Lyon (SLy4) equation of state (jChabanat et al.l 119981 : 
ISteiner fc Watts! l2009f) FI We find that Ncs/r is nearly 
constant over a large region (0 < r < 0.6R), which we 
expect is a feature of any NS equation of state (although 
the magnitude of Ncs/r will depend somewhat on the 
equation of state; see ? l3.3p . 

A NS can support various types of oscillation modes, 
including p- modes and (/-modes. The buoyancy that 
allows core g-modes to propagate is either thermally 

^ The model, including the crust, was kindly provided to us by 
A. Steiner. The p-mode and core g-mode that we consider couple 
in a region deep below the crust (see Fig. (3)1 ■ Our results are 
therefore not sensitive to the detailed properties of the crust and 
for simplicity we assume that the NS is a completely fluid body 
and that A''^ = in the crust (the outer ~ 1 km of the NS). 



WEINBERG, ARRAS, & BURKART 



0.1 



0.01 



10 



10" 



-3 







1 ' 

fs 








P 








- 








- 


Nc^r 




. 

\ •' ■ - 


: 
- 








. ^■ 
\ 1- 


- 

- 








: 1- 
1 h 












- / 










-1 










1 










\ 










I 














1 







10 



(km) 



Fig. 2. — Radial structure of an M = IAMq NS assuming the 
SLy4 equation of state: the hues show p (dash-triple-dot line; in 
units of lO'^^ g cm~^), Cs (dash-dot line; in units of ujqR where 
OJQ = (GM/R^y^^) , A'' (dash line; in units of ) , and Ncs/r (solid 
line; in units of i^g)- The dotted line shows that N/uiq ~ 0.06r/R 
for r < 5 km. 



induc ed (jMcDemiott et al.l 119881 : iGusakov fc Kantod 
l2012f) or due to pr oton-neutron composition gradients 
(jReisenegger fc Go ldrcich 1992). Here we focus on the 
latter and assume a normal fluid NS. Note, however, that 
at temperatures T < 10^ K, the bulk of the NS interior is 
expected to be a superfluid. The seismology of a super- 
fluid NS can differ significantly from those of a normal 
fluid NS. In particular, the core of a zero temperature 
superfluid NS does not ap pear to support pro pagating 

modes (Lee 1995; A ndersson fc Comeii l200l[) . How- 
ever, the core of a finite temp erature superfluid NS doe s 
support propagating (/-modes (jGusakov fc Kantoill2012[ ). 
albeit with properties that are sensitive to the equation 
of state, the core temperature, and the model of nucleon 
superfluidity. Given the uncertainties and difficulties as- 
sociated with modeling g-modes in a finite temperature 
superfluid core, we assume, for simplicity, a normal fluid 
N S whose core q-modes are simil ar to thos e described 
in iReisenegger fc GoldreichI ()1992[ ) and Eil (1994; and 
like these studies, we ignore general relativistic effects on 
the oscillations). Nonetheless, it is important to keep in 
mind that the results we present might be sensitive to 
superfluid effects. We briefly return to this topic in the 
conclusions (§ E]). 

In Figure[3]we show eigenfunctions of our NS model for 
a (p, g) daughter pair that satisfy equation ([1} in the core. 
(We normalize the eigenvectors a as oj^ j x p\a\^ — Eq, 
where p is the stellar density and Eq = GM^/R.) Be- 
cause their wavelengths match over such a large region 
(r < 0.6i?), the coupling between the modes is very 
strong. In Figure [3l we show the coupling of this pair 

to the equilibrium tide, k^^\ and in Figure 2] we show 
their coupling to an eigenmode parent, Kabc-, for parents 
that range from high-order (/-modes to the /-mode to 
low-order p-modes (we calculate the coupling coefficients 




r / R 

Fig. 3. — Radial profile of the eigenfunctions of a p-mode (solid 
grey line) and a g-mode (dashed line) for the SLy4 NS model; 
their degree, order, and frequency in units of ujq are (^ 5 , 7i„ , t*;^ ) = 
(1, 105, 200) and (£c, nc,ujc) = {i, -80, 6 X 10"*). The curves show 
bri^b and CftAcOJc, respectively, in units of Rloq. Note that the 
radial wavelengths of the modes are almost exactly equal for r < 

0.6R. The black line shows k^^''' (< r) for the coupling of the modes 
to each other and the equilibrium tide. 

using the expressions given in WAQB; see also §[211] be- 
low). Because the p-g coupling occurs well-below the 
crust (which resides at r > 0.9R), we do not expect it to 
be sensitive to the properties of the crust. 

Since Kabc is symmetric in abc, ii g : pg coupling can 
be strong then p : gg coupling can also be strong. For 
P'.gg coupling, the constraint \kb — kc\ ka implies that 
the daughter frequencies must satisfy AcWb ~ AfcWc (un- 
like with a (p, g) daughter pair, the daughter and parent 
wavelengths do not need to be equal) . Figure |4] shows 
Kabc for a self-coupled g-mode daughter coupled to a p- 
mode parent and demonstrates that p : gg coupling can 
indeed be strong. While there may be stellar systems 
where p : gg coupling is important, for the remainder of 
the paper we focus on parent waves that are either g- 
modes or a tidal perturbation, and therefore we do not 
further consider p : gg coupling. We also do not consider 
the coupling of the equilibrium tide to the dynamical tide 
and a p-mode even though it is a form of p-g coupling 
(since the dynamical tide is effectively a driven g-mode) . 
WAQB showed that such coupling can lead to a linear- 
like steady-state driving of the p-mode (see WAQB's sec- 
tion 9 and Appendix B.3); it is therefore very different 
from the nonlinear exponential driving that characterizes 
the PGI. Finally, for reasons given in § 13.11 we find that 
g'.pp and p.pg coupling are much weaker than the other 
types of p-g coupling. 

We now describe the properties of the p-g coupling 
coefficient in more detail. In § 13. II we present an analytic 
estimate of Kabc, in § 13.21 we apply this estimate to the 
NS model, and in § 13.31 we discuss the influence of the 
NS equation of state on the coupling strength. In § 13.41 
we explain why nonlinear terms beyond second-order in 
perturbation theory are not necessarily significant at the 
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PGI threshold. 

3.1. Properties of Kabc 

To see how Kabc depends on the stellar and mode pa- 
rameters, assume that the daughters are short wave- 
length, low degree modes. For a (6, c) = {p,g) daughter 
pair, this implies that in the propagation regions br 3> 6/i 
and Cr Ch, where the subscript r and h refer to the 
radial and horizontal displacement. 



J, \ f , CsSm(/)f, 

{br,bh)~ — coscpb, 

OJb \ ujbr 

, . Ac f LOcSUKpc COS0 
{Cr,Ch) ^ , — r— 



kb.cr is the rapidly varying phase, Ab^, 



(8) 
(9) 



= c7\Ic7'dr)- 



and 



[Eoab,c/pr^] ' , ab 

(A^/r)(/7Vdlnr)-i (see, e.g.. lAerts et~an[20Tol ). An or- 
dering of terms in WAQB's final expression for Kabc (lines 
A55-A62) reveals that their lines A60-A61 dominate and 
thus 



dKabc Pr'^ rp 2( N , 

-7j ^ -—FbUJ^ [Or - OhjbrCh 

dmr ho 

17 / ^l/2 

Pb [abac) ' , ^ 1^6 



(10) 



A, 



(Or — Oh) cos (\)b cos (/)c, (11) 



where Fb is an angular integral (see A21 in WAQB). An- 
gular momentum conservation requires that the modes 
satisfy the selection rules \lb — < < ^6 + with 
+ ^6 + even and rua -I- rub mc = 0. The angular 
factor Fb / Ac is typically of order unity and depends only 
weakly on the angular degree of the modes. 

In order to obtain a form for Kabc that allows for ac- 
curate numerical integration, WAQB performed a series 
of integration by parts on the origina l, compact form for 
i^abc (see also lWu fc Goldreichll200ll ). We took this nu- 
merically useful but non-compact form for Kabc as our 
starting point in deriving equation (|lip. However, for 
p-g coupling, it is straightforward to derive this same 
equation starting from the original, compact form. In 
Appendix [B] we show this for the particular case of the 
equilibrium tide coupled to a {p, g) daughter pair. 

The coupling scales as ujb/i^c ^ 1, the ratio of the 
p-mode to g-mode frequency. If the radial wavelengths 
match, LOb/ijJc — tc/tb, where tb^c — kb,crlu)b,c is the local 
radial group travel time of a mode. Physically, the larger 
the ratio uJb/oJc is, the longer the daughters (the (/-mode, 
in particular) spend in the strong interaction region and 
therefore the larger Kabc is. 

The condition that the wavelengths match (eq. [1]) im- 
plies that Kabc (X LUb/ujc c<: uj'^ (X iOc"^. In principle, Kabc 
can therefore be arbitrarily large. However, if we require 
that the waves are global normal modes (i.e., standing 
waves), there is an upper limit, uJb,ma,-x, to the p-mode 
frequency and a lower limit, a;c,ininj to the g-mode fre- 
quency. Depending on details of the stellar structure, 
W{,,max might be set by, e.g., the acoustic cutoff frequency 
of the atmosphere or the critical wavenumber kb oc LOb 
above which linear damping near the stellar surface is 
so rapid that the mode does not reflect. Similarly, since 



]we describe the physics that determines Wh^max and 
Wc,min for a NS. 

For g:pp coupling, Kabc oc i^b&r took b = c because 
the interaction is maximized for self-coupled daughters) 
while for g : pg coupling of equal wavelength daughters, 
Kabc oc Lu'^brCh (sec cq. [n])- Since br/ch ~ oOc/^b, g-PP 
coupling is weaker than g : pg coupling by the ratio of the 
g-mode to p-mode frequency; by symmetry, the same 
argument applies to p'.pg coupling. Because they are so 
much weaker, we do not further consider g:pp and p'-pg 
coupling. 

3.2. Neutron star core 

As FigureHshows, for r < 0.5i? of the SLy4 NS model, 
p w 10^5 g cm-3^ N/r w OMujq/R, and c, w 1.5i?a;o, 
where ojq = {GM/R^Y''^ ~ 1.07 x 10" rad s~^ is the 
dynamical frequency. Furthermore, we find ab ~ 0.6/ R 
and ac ~ OA/R. By equation ([1]), this implies that the p- 
mode and g-mode wavelengths are approximately equal 
if 

LJb — — 90Ac LjQ. (12) 

rujc \ ujc J 

This relation will inform our choice of ujb /i^c in the esti- 
mates below. 

For a high-order g-mode parent coupled to a {p,g) 
daughter pair with matched wavelengths, the coupling 
occurs near the parent's inner turning point ra — 
(aja/0.06a;o)i?, i.e., the location where Wq ~ N. This 
is because the product of all three waves is largest and 
nearly constant near Tq. From equations (|11|) and (|12p 
we find that at Tq, 

dKabc ^ O.OlFb / M UJb /c^o 

dlnr AaAc \pR^ J coc \i^a 

:^5xl0'(^^y(^^y\ (13) 
V200a;oy VlOOHzy ^ ' 

This expression agrees well with the magnitude and uja^ 
scaling of the full Kabc integration (see Figure H]) . 
If the parent is the equilibrium tide. 



GM 



gR \RJ 



r \ ' 



ah 



A2 



(14) 



where g is the gravity and x = £a — dlng/dlnr (see 
equations A12 and A13 in WAQB; although we include 
the gravitational perturbation due to the equilibrium tide 
in the full numerical calculations, we ignore it in the 
analytic estimates below because it is a small effect). In 
the NS core g ~ AirGpr/S and for £a = 2 and equal 
wavelength daughters we find from equations (jlip and 

m 



dK 



(oq) 
be 



din) 



3 X 10"* 



UJb 



V200wr 



R 



(15) 



The coupling scales linearly with r in the region where 
the daughter wavelengths match, in good agreement with 
the fuU K^^'^^ calculation (see Figure [3|)|3 

* WAQB showed that there is an additional, potentially impor- 



kc oc W-\ local damping might also determine Wc,min- In contribution to k^^^' from the linear inhomogeneous terms in 
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fa (Hz) 

Fig. 4. — Magnitude of the three mode couphng coefficient Kabc 
for a fixed daughter pair (fe, c) as a function of parent mode 
eigenfrequency fa = aJa/Svr. The parents are <?a = 2 modes 
and span the range from g-modes to p-modes (the /-mode is at 
fa = 2.5 X 10^ Hz). The filled circles show g .pg coupling assuming 
the daughters used in Figure \3\ The open circles show p : gg cou- 
pling for a self-coupled daughter that is the g-mode daughter used 
in Figure |3] The points extend slightly into the regimes o{ p: pg 
and g .gg coupling. The crosses show the inverse of the dynamical 

tide amplitude (see eq. [22]). 

The coupling strength is near the maximum values 
given by equations (IT^ and ([TSl) as long as |fcf, — fed ^ ka- 
Numerically, we find that high-order p-modes and g- 
modes satisfy the dispersion relations Wf, ~ iuj^rn, and 
LJc — 0.035cJoAc/nc. It follows that within the NS core, 
\kb — kc\ ^ 27r/i? for \nh — nc\ less than a few; thus, since 
the lengthscale of the equilibrium tide is « i?, for each p- 

mode, there are a few g-modes for which hi^^^^ is close to 
the maximum value. As we describe in § 21 one must ac- 
count for this effect when determining how many modes 
are unstable to the PGI. 

What is the maximum value of ljJjj/uJc for which equa- 
tions (fT5)) and (fTS)) are valid? In § 14.41 we show that the 
PGI requires that the p-modes (but not the ij-modes) 
form standing waves. We therefore argue that W(,,max 
is set by the acoustic cutoff frequency of the NS at- 
mosphere Wac — Cs/2Hp, where Hp is the density scale 
height. For w;, > Wac, acoustic waves do not reflect at 
the stellar surface and form standing waves. In Ap- 
pendix [C] we show that for a cold NS, cjac — 2OO/94 " '^cjo, 
where p4 = p/lO"* g cm^^. At temperatures T ~ 10^ K 
and densities p < 10^ g cm"'^, ideal gas pressure dom- 
inates and Wac ~ IOOOwq. We therefore expect the NS 
acoustic cutoff frequency to lie somewhere in the range 

100 < Wac/t^O < 1000. 

Linear damping can also potentially limit the maxi- 
mum value of ujb/ujc- In particular, if the p-mode lin- 
ear damping rate exceeds the reciprocal of its round-trip 
travel time between turning points, ttravei — 2 / dr/cg ~ 
2/ti;o , it will decay before reflecting and forming a stand- 

thc equations of motion (see their Appendix A). For a (p,g) daugh- 
ter pair, however, the homogeneous terms dominate the coupling. 



ing wave (see, e.g.. iGoodman fc Dicksonlll998D . In Ap- 
pendix [C] we calculate the linear damping rate of modes 
and show that for Uf, < Wac, the linear damping rate 
7f, <C ^t^avoi- (Although the PGI does not require the g- 
mode to form a standing wave, we note that jc < ^travel 
for low-degree (7- modes with frequency uJc ^ 10~^aJo: 
where here ^travel — '^-i^c/^^t) J Ndlnr). Moreover, for 
low-degree modes that are unstable to the PGI (such as 
those shown in Figures |3] and |4]) , we find jb.c ^ ^b,c and 
thus the modes are well-described by the adiabatic stel- 
lar oscillation equations. We therefore conclude that for 
a (p, g) daughter pair with nearly equal wavelengths (i.e., 
satisfying eq. 1121). the maximum value of uj^/lOc is 

„,ax(2)«4xlO=A-(5^)'^ (16) 

This result motivates taking — 200a;o as a reference 
value in equations (IT51) and (IT5|) . 

3.3. Influence of the equation of state on Kabc 

We base our calculations of Kabc on a single type of nu- 
clear equation of state (SLy4). While SLy 4 is consistent 
with all current observational constraints (jSteiner et al.l 
our choice is otherwise arbitrary; ideally we would 
like to compute Kabc for different equations of state. 
Unfortunately, most microscopic calculations of high- 
density nuclear matter do not provide sufficient details 
to enable calculation of the buoyancy iV c>c — Cg and 
therefore Kabc- In particular, they often provide the equi- 
librium sound speed Ce but not the adiabatic sound speed 
Cs , where adiabatic in this context implies constant com- 
position fReisenegger & Goldreich 1992; Lai 1994). 

Although we do not have a precise estimate of N 
for other equations of state, w e believe that the results 
for SL y4 are representative. iReisenegger fc GoldreichI 
(|1992D showed that in a NS core N « (xjifl'^glcc, where 
X is the proton fraction. Because p is nearly constant over 
a large fraction of the core , x^^'^ and Ce — Cs change very 
little with r while g (x r (between the stellar center and 
r ~ ^/2). Therefore, for any equation of state, we ex- 
pect Ncs/r to be nearly constant over a large portion of 
the star, i.e., similar to what is shown in Figure [2] This 
implies that the p-mode and g-mode wavelengths can 
be equal over a large region, allowing for a strong non- 
li near coupling. Furthermore, ba sed on the discussion 
in IReisenegger fc GoldreichI (|1992D and the SLy4 results, 
we expect the overall magnitude of Ncg/r to be within 
a factor of a few of the SLy4 value for any viable equa- 
tion of state0 Because Ncs/r determines the ratio uJb/oJc 
and thus Kabc (see eqs. [TT] and [H]), we conclude that 
most equations of state likely yield qualitatively similar 
p-g coupling results. 

3.4. Higher-order p-g coupling 

^ ILail 111 9941 ) apply an approximate fitting method to extract 
N fro m the four micro scopic equations of state of Wirin ga et al.l 
{1983). Figure 3 in ILail shows that all four yield high-order g-mode 
frequencies that are within a factor of two of each other for a given 
radial order. Moreover, three of the four have g-mode frequencies 
almost exactly equal to that of th e SLy4 model (cf., the g-mode 
dispersion relation given in ^ 13.21 1. Since uJc <x A^, this suggests 
that N is likely within a factor of fx 2 of the same value in all five 
equations of state. 
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As noted in § [21 at the PGI threshold the hnear terms 
are comparable in magnitude to the lowest order (n = 2) 
nonlinear terms. Since the standard ordering of terms 
does not apply when going from n = 1 to n = 2, one 
might wonder whether it is valid to neglect n > 2 terms. 
By carrying out a rough estimate of the magnitude of 
the n = 3 terms, we now show that higher-order terms 
are unlikely to be important near threshold. 

At order n = 3, the amplitude equation of a daughter 
mode b includes four-wave coupling terms of the form 
KMefqdIelfj where Kbdef is the four- wave coupling co- 
efficient. The amplitude equations of the other modes 
include analogous four- wave terms. There are four types 
of four-wave couplings to consider: (i) a ^ {d,e,f}, (ii) 
a = d ^ {e, /}, (iii) a = d = e ^ /, (iv) a = d = e = f. 
The last two cases correspond to a self-coupled parent. 
Because we are interested in assessing whether higher- 
order terms influence the onset of the instability, we as- 
sume the parent is at its initial, linear amplitude and the 
daughter amplitudes are infinitesimal (as in the stability 
analysis of § [2]) . 

In cases (i) and (ii), the n = 3 term in the ampli- 
tude equation contains the product of three daughter 
amplitudes and two daughter amplitudes, respectively, 
whereas the n — 2 term contains only a single daughter 
amplitude. Since the daughter amplitudes are infinites- 
imal, the case (i) and (ii) n — 3 terms will necessarily 
be negligible compared to the n = 2 terms at threshold; 
therefore they cannot prevent the onset of the p-g insta- 
bility (they might influence the saturation, however; see 

In case (iii), the n = 3 term contains the product 
of two parent amplitudes q^, whereas the n = 2 term 
contains a single parent amplitude qa- However, the 
argument used in cases (i) and (ii) does not immedi- 
ately apply because the parent amplitude is finite at 
threshold. Instead, to determine whether n = 3 terms 
are important near threshold, we need to determine if 
\qaKaabf\ > \Kabc\- Notc that the forn i of Knhr.d is simila r 
to that of Kabc (both are derived in I Van HoolstI 1199^ : 
in particular, the integrand of Kabc contains terms of 
the form (V • l)^, V • ^S^j^^,, and C;jCffeCfi while the in- 
tegrand of Kabcd contains terms of the form (V • ^)^, 
(V • i)^^,e,, (e^,C:^)^ and ^.ek^^^^^- Here we use C 
to represent the Lagrangian displacement of each of the 
modes (e.g., (V • = V • aV • 6V • c), the subscript 
semi-colon denotes a covariant derivative, and we did not 
write down the terms involving the gravitational poten- 
tial and its perturbations because they are negligible for 
the coupling of high-order modes. Comparing forms, we 
see that the terms in Kaabf contain an extra factor of 
the parent's spatial derivative ~ dai/dx^ relative to the 
terms in Kabc- Assuming kf c^i kb (otherwise Kaabf is neg- 
ligible for the reasons given in § [Ij , then since kb — kc we 
see that the n = 3 terms are important near threshold if 
\qadai/ dx^ \ > 1 (note that this statement is independent 
of our choice of normalization) . 

For a parent that is a low-order mode or the equi- 
librium tide, dtti/dx^ < 1. Since we are interested 
in cases where \qa\ ^ 1 (e.g., for the equilibrium tide 
\qa\ {M'/M){R/af < 1), we have \qada,/dxi\ < 1 
and thus the n — 3 term is negligible compared to the 



n — 2 term near threshold. 

For a parent that is a high-order 5-mode (e.g., a dy- 
namical tide mode), dai/dx^ ^ kaUr 3> 1, where ka is 
the radial wavenumber of the parent. Thus, n = 3 terms 
are comparable to n — 2 terms if IqafcaOrl ^ 1- This is 
just the usual nonlinearity parameter that determines the 
critical threshold above whic h a q-mode overturns the lo- 
cal stratific ation and breaks ([Goodman fc Dicksonlll998l : 
[Barker fc~b gilvic 20l3). If Kabc > max(fcaar), there is a 
range of |qa| where p-g coupling is unstable (|«;abc(Za| > 1) 
but the n = 3 term is negligible {\qakaa.r\ < !)■ Us- 
ing the values for the NS core and equation ([U, we find 
msxikadr) ~ 4 X 10~'^Aa(a;o/a;a)^, with the maximum 
occuring at the mode's inner turning point. Combining 
this with our estimate of Kabc for g'-pg coupling (eq. |13| : 
see also Fig. [IJ, we find that for uja/^o ^ 10~^Aa, there 
is a range of \qa\ for which the parent is p-g unstable 
but the n = 3 terms are negligible. For the particu- 
lar case of the dynamical tide in coalescing NS binaries, 
\qakaar\ ^ 1 at all frequencies fa (see eq. [22] below) 
and therefore, as is the case for the equilibrium tide, the 
n = 3 terms are negligible at the PGI threshold. 

The analysis of case (iv) is similar to that of case (iii) 
except now the n = 3 terms are important near threshold 
if \qadai/dx^\^ > 1. Since we found \qadai/dx^\ <C 1 in 
analyzing case (iii), the n = 3 terms of case (iv) are also 
negligible near threshold. 

Similar arguments apply to yet higher-order terms. We 
therefore conclude that even though the magnitude of the 
n = 1 and n = 2 terms are similar at the p-g instability 
threshold, the n > 2 terms are not necessarily significant 
and a perturbative approach remains valid. 

4. INSTABILITY IN COALESCING NEUTRON 
STAR BINARIES 

In this section we consider the stability of the tide in 
coalescing NS-NS and NS-BH binaries. In § § [O and [U 
we determine when the equilibrium tide and dynamical 
tide are unstable to the PGI, respectively. We then com- 
pute the nonlinear growth rates of the unstable daugh- 
ters. For comparison with the PGI, we also consider the 
stability of the tide to the resonant parametric instabil- 
ity. In § 14.31 we consider the implications if the parent is 
well-above the PGI threshold. In 14.41 we show that the 
PGI growth rates are so large that the g-modes (but not 
the p-modes) grow significantly in less than their group 
travel time across the star, implying that the g-mode 
driving is local. 

4.1. Stability of the equilibrium tide 

We derive the equilibrium tide p-g stability criterion 
in Appendix [^ For a circular orbit, the dominant £ = 2 

equilibrium tide is unstable to the PGI if e ^ I'^ic'*'' 1""'^ 
(cf. eq. [5]), where the tidal amplitude factor 

\MtJ viooHzy ^ ' 

Here M' is the companion mass, Mt = M-\-M' is the to- 
tal mass of the binary, and /gw is the gravitational wave 
frequency (which equals twice the orbital frequency). 
From equation ([T5|) . we thus find that the equilibrium 
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tide is unstable to the PGI when 

When unstable, the equilibrium tide excites daughter 
pairs that grow exponentially at a rate 



(18) 



:75 



M'\ f 



UJb 



MtJ V 200^0 



gw 



100 Hz 



Hz. 



(19) 



This is much larger than the linear damping rate of £ = 2 
modes in the frequency range of interest (see Appendix 
[C|) . Moreover, it is much larger than the inverse of the 
gravitational wave inspiral timescale 



"-gw - 



/gw 
/gw 



5.9 



M 



-5/3 



/g 



100 Hz 



-8/3 



s (20) 



_1.2M0, 

(jPeters fc Mathewill963f ). where the chirp mass M. = 
[{MM')yMt]^/^ (for M = M' = IAMq, M ~ 1.2Mo). 
The number of e-foldings that the daughters can grow 
before the binary merges is 



ao^ 



^(eq) c(/gw 
he J. 

/gw 



2. OA/, 



o 



25 Hz 



,(21) 



where fi is the value of /gw when the daughters first be- 
come unstable as given by equation (jl8p . Since the insta- 
bility does not rely on resonant interactions, the daugh- 
ters are continuously driven as the binary inspirals. The 
scaling /"^/'^ implies that most of the growth occurs at 
low frequencies, where the orbital decay is slowest. Be- 
cause the number of e-foldings is so large, even daughters 
with very small initial amplitude can reach a significant 
amplitude well before the binary merges; the maximum 
amplitudes are therefore set by the nonlinear saturation 
of the instability rather than the time until merger. 

The preceding calculation is for a single daughter pair. 
However, since rib — 100 for uJb — 200cjo (see the p-mode 
dispersion relation given in § 13. 2p . there are « 10 dis- 
tinct p-modes that couple to the equilibrium tide with 
near equal effectiveness. Furthermore, as noted in § 13. 2[ 

for each p-mode there are a few (7-modes for which n^'^^^ 
is near the maximum value. This suggests that for a 
given {£b,lc), there are ss 10 — 100 daughter pairs that 
have similarly large values of kI^^\ And since the mag- 
nitude of is a weak function of the daughters' angu- 
lar degree (it decreases only slightly with £b,c), the total 
number can be greater still. The number of daughters 
pairs that are p-g unstable to the equilibrium tide can 
thus be Ncq > 100. Moreover, the number increases as 

the orbit shrinks and pairs with ever smaller be- 
come unstable. The net rate at which the PGI dissipates 
the equilibrium tide's energy might therefore be consid- 
erably larger than the rate due to a single daughter pair 
(see SIOI). 

We apply a similar analysis to evaluate the stability 
of the equilibrium tide to parametric resonance. This 



involves the coupling of the equilibrium tide to a pair 
of g-mode daughters with small detuning relative to the 
tidal frequency, |A| = \uJb + ujc - wtidd < Wtido (for a 
circular orbit Wtide equals twice the orbital frequency). 
Such daughters are unstable if their nonlinear growth 

rate Tbc > a/A^ + jbjc, where Tbc ~ Wtiddetic'^^ | (see 
WAQB). We find that for self-coupled 5-mode daughters 



,(<=q) 



0.3, nearly independent of mode frequency and 
angular degree^ The parametric growth rate is therefore 
Tbc ~ 10~'^(/gw/10 Hz)'^ Hz for an equal mass binary. 
While this is larger than 7 for resonant low £ modes with 
frequencies > 10 Hz (see Appendix . it is smaller than 
their average detuningQ We therefore conclude that the 
equilibrium tide is stable to parametric resonance. 



4.2. Stability of the dynamical tide 
The dynamical tide is unstable to the PGI if \Kabc\ 



> 



\qa (cf. eq. [5]), where |ga'^"^"''| is the amplitude 

of the dynamical tide mode after it has under gone linear 
resonan t driving. Following the calcula tion bv lLail (1994; 
see also IReisenegger fc Goldreichlll994[ ). we find 



3 X IQ-'-'h 



fa 



100 Hz 



19/6 



(22) 



where h = {W /Mf/'^{2M /Mtf/^ and we made use of 
the expression for the linear overlap integral given in Ap- 
pendix [C] (note that o ur eigenfu nction normalization is 
different than that of iLail [l994l ). The mode frequency 
fa = L^al^T^ cquals the gravitational wave frequency /gw 
when the resonance occurs. We verified this estimate of 
|^(dyn)| numerically integrating the linear amplitude 
equations for an orbit decaying due to gravitational ra- 
diation. 

Since the dynamical tide mode a is a g-mode, it couples 
to a (6, c) = (p, g) daughter pair with a strength Kabc 
given approximately by equation ()13p . We thus find that 
the dynamical tide is unstable to the PGI if 



/a > 12 h 



7 / '^h 



V 200^0 



-12/7 



Hz, 



(23) 



in good agreement with the full numerical calculation 
shown in Figure |4l The result depends only weakly on 
the mass ratio M'/M; e.g., for a NS-BH system with 
M = I.4M0 and M' = IOM0, the dynamical tide is 
unstable for fa ^ 14 Hz. 

When unstable, the dynamical tide mode excites 
daughter pairs that grow exponentially at a rate 

l^ic '^^"'' — '^c|Ka6c9a'^^"''|, whcrc the superscript (a,dyn) 



,(<=q) 



The arguments given in section 5.3 of WAQB explain why 
_ 0.3. The coupling is weaker than that of daughter 
pairs by the factor maxjixiij/tJc} ~ 10^ (see eq. 1161 ). Note that we 
focus on self-coupled pairs because, as WAQB showed, the coupling 
strength peaks strongly for self-coupling. 

^ The g-mode dispersion relation given in § 13.21 determines the 
frequency spacing of the modes and implies that for self-coupling, 
A| ~ (10/Aa)(/a/10 Hz)^ rad s"-"-. This is always greater than 
the resonant parametric growth rate V^y^ for low la modes. As 
the binary inspirals and a;tide increases, |A| will be smaller than 
average for brief intervals. However, because V\y^ is so small, the 
daughters will not have a chance to grow significantly during these 
intervals. 
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indicates that mode a is a dynamical tide mode that was, 
at some earUer time, resonantly excited by the tide (see 
eq. [6]). From equations (|13l) and (|22|) . we find 



^(a.dyn) 
be 



72h 



V200wo 



fa 



100 Hz 



7/6 



Hz. 



(24) 



Like the daughter growth rate r["''^'^"-* is much 

larger than and the linear damping rate of relevant 
i = 2 modes. The number of e-foldings that the daugh- 
ters can grow between the time when they first become 
unstable at /gw = fi and the binary merges is 



M 



l.2Mp 



-5/3 



20 Hz 



-3/2 



(25) 



where we used a value for fi motivated by equation (1^^ . 
This estimate assumes that nonlinear interactions dp not 
prevent the parent from reaching the amplitude \qa | 
in the first place. Whether this is true depends on the 
initial amplitude of the daughters and the duration of the 
parent's linear excitation; if it is not true, the growth rate 
and the number of e-foldings will be smaller. 

As in the case of the equilibrium tide, there is enough 
time for many e-foldings of growth before the binary 
merges, with most of the growth occurring at low fre- 
quencies. And like the equilibrium tide, there can be 
many daughter pairs that are unstable to the dynami- 
cal tide, i.e., A'^dyn ^1- A key difference, however, is 
that the equilibrium tide is continuously driven as the 
binary inspirals whereas the dynamical tide modes are 
only driven significantly for a brief interval near their 
linear resonance. We discuss the consequences of this in 

m 

We now evaluate the stability of the dynamical 
tide to parametric resonance. The maximum growth 
rate of parametrically unstable daughters is Fbc ~ 

Wa|Kabc9a"^'*'"'' I, where Kabc IS the coupling coefficient for 
the dynamical tide mode a coupled to pair of g-mode 
daughters (6, c) with frequencies Wf, -I- Wc — Wa . The rel- 
evant parent frequency is Wa rather than wtidc because 
post-resonance, the dynamical tide mode oscillates at its 
natural frequency uja (^Lai 1994). For fa ^ 10 Hz, we 
find Kabc ~ 100 nearly independent of fa- The growth 
rate is thus Tbc ~ 5h{fa/ 100 Hz)^^/^ Hz, which is larger 
than the linear damping rate of low £ daughters with fre- 
quency > 10 Hz (see Appendix [C]) and their minimum 
A I (because the dynamical tide Kabc is large as long as 
iT-b — ricl ria, 1^1 '^^n. be much smaller than the self- 
coupled |A| used in ? l4.2l to determine the stability of the 
equilibrium tide to parametric resonance). We therefore 
conclude that the dynamical tide is unstable to paramet- 
ric resonance for fa > 10 HzEl However, since the PGI 
has a much larger growth rate for 10 < fa ;$ 200 Hz, it 
is more likely the dominant instability. 

4.3. Validity of the PGI analysis well above threshold 



® For fa < 10 Hz we find K^bc ~ 100(/a/10 Hz)-^ (it increases 
as because at low frequencies ar oc r^'^ oc fa'^, where Va is 
the inner turning point). Nonetheless, the dynamical tide is stable 

for fa $ 10 Hz due to the strong frequency dependence of Igi"^^"'! 
and the linear damping rate 7. 



Near the PGI threshold, the linear forces that act on a 
perturbation are comparable in magnitude to the n = 2 
forces (and much larger than the n > 2 forces; see §[2]and 
S 13. 4|) . However, well-above threshold the n = 2 forces 
dominate and it is not clear if the linear relations that 
we use to analyze the p-g coupling remain valid. In par- 
ticular, if the restoring force acting on a perturbation is 
significantly affected by nonlinear interactions, is it valid 
to use, e.g., the linear eigenfunctions and dispersion rela- 
tion, to evaluate quantities like Habc and the wavelength 
matching criterion fcf, ~ fee? Mathematically, one is al- 
ways free to choose the linear basis in order to expand 
the perturbations, so perhaps this is not a concern (at 
least during the early stages of the PGI when the n > 2 
forces remain smaller than the n = 2 forces and the per- 
turbation series converges). On the other hand, it is not 
clear what it means to label a perturbation a p-mode 
or a g-mode when the linear restoring forces that give 
these modes their properties are overwhelmed by nonlin- 
ear forces. 

Although this is an important concern for the PGI in 
general, it is unlikely to significantly impact our main 
conclusions regarding the infiuence of the PGI on coalesc- 
ing NS binaries. As we describe in § [SJ the orbital phase 
error and tidal heating due to the PGI occur primar- 
ily near the equilibrium tide instability threshold (i.e., 
near fi). The linear relations that we use to analyze the 
PGI should therefore be approximately correct during 
this most important phase of the PGI. 

4.4. Local driving 

The PGI growth rates (eqs. [19] and [24]) assume that 
the daughters are global standing waves. If, however, a 
daughter's growth rate within some region is much larger 
than its inverse group travel time across that region, the 
daughter undergoes runaway local growth. It is then 
more proper to treat the daughter as a traveling wave 
rather than a standing wave. Although this does not 
affect the stability criterion (near threshold the growth 
rates are necessarily small), as we describe in § (5] it in- 
fluences the saturation of the instability. 

In Appendix |D] we show that the local growth rate of 
the PGI is 



fbcir)^uJc[abac] 



dHa 



be 



dr 



Qaittr - ah) 



(26) 
(27) 



where we use the symbol F to indicate a local growth 
rate and the second line follows from equation (fTTj) . 
Based on the results of § 13. 2[ we see that for the equilib- 
rium tide the local growth rate is larger than the global 
growth rate by a factor of about {R^abUc)"^ ~ 2, nearly 
independent of radius (recall that for the equilibrium 
tide Qa — > e). The local and global growth rates are 
similar because the coupling is approximately constant 
throughout the core, which extends out to r « R/2. 
For the dynamical tide, the local rate is larger than 
the global rate by « {r/R)~^ because the coupling is 
strongest at small radii where the dynamical tide peaks; 
the maximum occurs at the tide's inner turning point 
r,/i? «0.1(/a/10 Hz). 
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For a p-mode, the radial group travel time in the core 
is tb = r/cs ~ 60{r/R)^s and we find 

Pr'-0.03(^)(^)5 (28) 

«0.01i(^^y'° (29) 

where for the dynamical tide we evaluated the result at 
Va- Since Tbctt <C 1 for both the equilibrium and dy- 
namical tides, the p-mode driving is in the standing wave 
limit. By contrast, for a g-mode with wavelength equal 
to that of a p-mode, tc ~ if,(iujf,/wc) ^ h (see § 13. ip and 
we find 

r<:«,.3xio3(^)(^)'l (», 

The y-mode driving is therefore well into the traveling 
wave limit for typical parameter values. This means that 
the g-mode grows significantly in the time it takes to 
propagate across a small fraction of the star (unlike the 
p-mode, which undergoes many reflections in that same 
time). 

5. SATURATION OF THE p-g INSTABILITY 

Having evaluated the stability of daughters to the PGI 
in the previous sections, we now consider their satura- 
tion. As we discuss in § 15.11 determining precisely how, 
and at what energy, the daughters saturate is a challeng- 
ing calculation that is beyond the scope of this paper. 
Instead, we leave the saturation energy as a free param- 
eter whose magnitude we attempt to constrain based on 
a few informed assumptions. 

We then use this result to determine how the PGI 
might influence the orbital evolution and heating rate of 
coalescing NS binaries. Because tidal interactions trans- 
fer energy and angular momentum from the binary orbit 
to the NS, they modify the rate of inspiral. This induces 
an orbital phase error A(j) relative to that of two point 
masses that accumulates over the course of the inspiral. 
Furthermore, the waves excited by the tide heat the NS 
interior and deposit angular momentum that can spin-up 
the star. In linear tidal theory, these effects are deter- 
mined by the amplitude and linear damping rate of the 
pare nts, i.e., o f the equilibrium and dynamical tides (see, 
e.g.. lLailll994[) . In nonlinear tidal theory, however, the 
daughters provide the parents with an additional, am- 
plitude dependent, source of dissipation. In § 15. 2[ we 
evaluate the nonlinear energy dissipation rate due to the 
PGI and from that result, estimate A(p in § 15. 3[ and the 
tidal heating and spin-up of the NS in § 15.41 

5.1. Daughter saturation energy 

The daughters' growth saturates when their nonlin- 
ear driving rate balances their nonlinear damping rate. 
The latter depends on whether the saturation occurs 
in the "discrete" limit or the "continuum" limit (see 
I Arras et"aI]|2003D . In the discrete limit, the daughters 
saturate by transferring their energy to a set of dis- 
crete modes whose linear damping rates are greater than 




6 8 10 12 14 



logio P (g cm ^) 

Fig. 5. — Nonlinearity parameter \ki,br\ and energy flux Ff, for 
the p-mode used in Figure [3] as a function of stellar density. The 
flux is in units of 10'^i?o'^o^~^- 



the daughters' nonlinear driving rate. In the continuum 
limit, the daughters saturate by exciting a turbulent cas- 
cade in which the energy input at the outer scale (that 
of the parents) cascades down, via nonlinear interactions, 
to a very short wavelength inner scale where it rapidly 
thermalizes. 

For coalescing NS binaries, the PGI growth rates are 
much larger than the linear damping rates of even the 
most highly damped, quasi-adiabatic g-modes (?? l4.1l and 
14. 2p . Moreover, since the driving of the daughter modes 
is local (§ 14. 4p . their saturation must also occur locally. 
Together these suggest that for this problem, the satura- 
tion occurs in the continuum limit and involves the de- 
velopment of a turbulent cascade. Because of the rapid 
local growth, the cascade might not saturate until the 
daughters become highly nonlinear and approach their 
local "wave breaking" amplitude |(?break(?')|l_f The daugh- 
ter amplitudes probably cannot signiflcantly exceed the 

linear amplitude of the parent Igi''"'*!, however. These 
considerations motivate us to parametrize the local sat- 
uration energy as 

i?,at(r) =(3mm (igbreakP, ) Eo, (32) 

where the magnitude of the parameter /? is unknown but 
may be of order unity according to the above arguments. 

Both modes and p-modes break when Ig/cr^^l 
1; at this amplitude, g-modes begin to overturn 
the stratification (see, e.g., Goodman & Dickson .19981: 
[Barker fc Ogilvij|2010f) and p-modcs impart order unity 
density perturbations. From equation we find that 
g-modes whose wavelengths match their p-mode pair 

" It is not clear whether |gbrcak(^')l ^ strict upper bound to the 
daughter amplitudes. Since the daughters oscillate at a frequency 
well-below their natural frequencies (ai, <g uji, and ac u}c', see 
Appendix , they might be able to reach amplitudes in excess of 

kbroakC?")!- 
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(fcc ~ kh) break at an amplitude 



|9c,brcak('')| ^ (kcCrY 



10' 



V 200^0/ V0.5i?/ 



(33) 

For p-modes, khbr peaks near the stellar surface and for 
ujb w lOOwo we find {hbr)'^ ~ 10"'' - 10"^ near the 
outer turning point (see Figure [5]). While this suggests 
that I gb, break I ^ I gc, break I, Figure [S]shows that only a very 
small fraction of the p-mode energy flux Ff, = uj^pr^b^Cs 
reaches the stellar surface (for a standing wave, its more 
proper to regard -Ff, as a measure of the mode energy 
density rather than flux). As a result, wave breaking 
near the surface might not saturate the p-mode growth 
occurring deeper within the star. Instead, the p-mode 
might saturate only when it breaks in the core, which 
from equation 1^ occurs at an amplitude (see also Figure 



kb,broak(?')| ^ [khbrY 



0.1 



0.5i? 



(34) 



In that case, since |(7c,broak| ^ |9b,broak| at all radii, the 
g-mode always breaks before the p-mode (recall from § [5] 
that \qb\ — \qc\ during the unstable growth). We will 
therefore assume that |gbroak| = |gc,brcak|- 

The amplitude of the equilibrium tide e > |gc, break | at 
all radii r < R/2 when /gw > 30{M' /Mt)-^^^ Hz (see 
eqs. [17] and [SS]). There is therefore enough energy in 
the equilibrium tide at these frequencies that the g-mode 
daughters can, in principle, undergo local wave breaking 
throughout the core (i.e., reach /3 ~ 1 at all r < i?/2; 
recall that the equilibrium tide driving rate is nearly in- 
dependent of radius). Since |(7c,break| r^, the outer radii 
r ~ i?/2 would then be the principal seat of energy dissi- 
pation. Based on equation p2[) . we therefore express the 
volume integrated saturation energy of daughters driven 
by the equilibrium tide as 



(cq) 



= /3cq|9c,brcak(^/2)| Eq 



(35) 
(36) 



where /3eq ~ 1 might be a reasonable estimate. 

For the dynamical tide, the driving occurs predom- 
inantly at the tide's inner turning radius ra/R ~ 
0.1 (/a/10 Hz). Plugging this radius into equation 
and comparing the result with equation (P^ . we find 

lo'a'^'^"^! ^ I ?c, break I at all Orbital frequencies; therefore, 
daughters driven by the dynamical tide do not break. In- 



< 



l^abc 



stead, since the instability turns off once |<7a| 
the daughters more likely saturate at an amplitude 



l^abc 



(37) 



i.e., after its resonant excitation, the dynamical tide 
mode a transfers nearly all of its energy to the daughters. 
Based on equation ([32]) , we therefore express the volume 
integrated saturation energy of daughters driven by the 
dynamical tide as 



E. 



(a,dyn) _ 



10" Vdyn/i' 



fa 



100 Hz 



19/3 



En 



(38) 
(39) 



where /3dyn ^ 1 might be a reasonable estimate. 

5.2. Nonlinear energy dissipation rate 
The nonlinear energy dissipation rate due to a single, 
saturated, daughter pair is Etc ~ 2rbci?sat (a factor of 
two appears because Pbc is the growth rate of the am- 
plitude not the energy). From equations ([19]) and ([35]) . 
the total dissipation rate due to the instability of the 
equilibrium tide is therefore 



E, 



eq^^lO^«/3eqA^eq(^ 



M'\ ( \ ' 

V200a;oy 



f k 



V 100 Hz 



erg s 



(40) 



where iVoq is the number of unstable daughters driven by 
the equilibrium tide (see ii l4.1l) and we took a single value 
of ujb as representative of all the unstable pairs. Similarly, 
from equations ([24]) and ([38]) . the total dissipation rate 
due to the instability of the dynamical tide is 



Edyn ~ 10^''/3dynA^dyn/i^ 



LOb 

200a;o 



fa 



100 Hz 



45/6 



(41) 



For a NS-NS binary at /gw — 100 Hz, these dissipation 
rates are ~ \(f N^qficq and ~ 10^A'^dyn/3dyn times larger 
than the linear dissipation rates of the equ ilibrium and 
dynamical tides, respectively fsee lLai|[l994t this assumes 
a NS core temperature T — 10^ K). Physically, this is 
possible because the PGI dissipates the energy in the 
tide on the rapid nonlinear driving timescale rather than 
the slow linear damping timescale. 

Since each dynamical tide mode is excited for only a 
brief interval during its linear resonance, the mode ceases 
to heat the NS (i.e., £;(°''^>'") ^ 0) once the nonlinear in- 
teractions dissipate all of its energy. By contrast, since 
the fluid always seeks to follow gravitational equipoten- 
tials, there is a continuous flow of energy and angular 
momentum from the orbit into the equilibrium tide. Un- 
stable daughters driven by the equilibrium tide therefore 
continuously dissipate energy at some finite rate E- 



^(cq) 

fee 



5.3. Phase error A(/) 

We now calculate the orbital phase error due to 
nonlinear tidal interactions. The total energy of the bi- 
nary system is i^tot = ^'orb + E, where i?orb is the orbital 
energy and E is the energy in the modes and the tidal in - 
teraction. Following a similar calculation by iLail (|1994D , 
the total power is 



^ dEtot J. I dE^rh 

-tjtnt = — — / = 



dE 



df \ df ' df 

and the change in phase over a time dt is 



/ 



(42) 



= 2TTfdt = 2nf^ 



df (dE, 



E,. 



^orb 



df 



The first term in equation (|43p is the point-mass result. 
The second terms represents the phase error due to the 
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tidally induced change in stellar energy 

,x ^ n df dE , df dE E ,^ 



where the gravitational wave luminosity 



E, 



(44) 



Fj — 1 



X 10' 



,51 



7W 



1.2M, 



10/3 



( /s 



V 100 Hz 



10/3 



erg s 



(45) 

and the second expression in equation (|44p follows be- 
cause -Etot — ^'gw 3> i?, i.e., gravitational wave emission 
is the dominant source of orbital decay. 

For the equilibrium tide, E is given by equation (|40p 
and we find upon integrating equation (|44l) 



-1 / r \ -3 



xf^yV^^y (46) 

(47) 

where in the second expression we used equation (jl8p to 
determine /i, the value of /gw at which daughters first 
become unstable. In § 14.11 we argued that N^^ > 100 
and possibly significantly larger. Thus, even if /3oq ^ 
10~^ — 10~^ (i.e. considerably smaller than the upper 
bound set by energetic arguments), the instability of the 
equilibrium tide can induce a phase error A(/>cq > 1. 
Moreover, since A0cq oc f~^, the phase error accumu- 
lates primarily at large orbital separations. 

As noted in § 15.21 nonlinear interactions do not affect 
the amount of orbital energy put into the dynamical tide 
mode, they only alter the rate at which the dynamical 
tide mode energy is thermalized (assuming again that 
the interactions do not influence the resonant excitation 
itself; see § 14. 2p . As a result, even if the dynamical tide 
is unstable to the PGI, its phase e rror is the same as that 
of linear theory, i.e., A^dyn < 1 (|Lailll994l ). 

5.4. Tidal heating and spin-up of the neutron star 

The nonlinear dissipation rates Ecq and i?dyn (eqs. [40] 
and [41]) describe how quickly orbital energy is con- 
verted to mode energy at saturation. Assuming that 
the saturation is steady (at least in a time-averaged 
sense), the modes deposit their energy into the star at 
a rate E — E^q + E^yn- Some of this energy will heat 
the NS and some of it will go into rotational energy. 
iBildsten fc Cutleil (|1992[ ) showed that the equilibrium 
tide cannot synchronize the spin of a standard NS. Al- 
though their calculation accounted only for the linear 
terms in the interaction Hamiltonian (for the nonlin- 
ear terms, see, e.g., WAQB), it is very unlikely that 
the PGI will synchronize the spin. This is because for 
spin frequencies /spin ^ 20 Hz the rotational energy 
Erot — 10^-^ (/spin/100 Hz)^ erg is larger than the to- 
tal energy removed from the orbit AE ~ J Eeqdt ~ 



10'*^(/,/100 Hz)-2/3 erg, where we assumed uniform ro- 
tation and used equation pO|) . From equation ([T8| . 
fi > 20 Hz and thus there is at most a small interval 
near the onset of the instability when the star can, in 
principle, briefly synchronize. 

While we conclude that synchronization is unlikely, 
without a detailed understanding of the saturation we 
do not know what fraction of E goes into heating the NS 
as opposed to spinning it up. For simplicity, assume that 
all of E goes into heating the NS. The thermal evolution 
of the NS heat content U is then 



— ^cq I -C^dyn i -'-^cool 



where 



C/~4.5 X 10"^ ' — ^ 



VlOS K 



erg 



(48) 
(49) 



(lLail[T99l and we neglect the cooling term i?cooi due 
to neutrino emission and surface phot on emiss ion which 
are sm all at temperatures T < lO^'^ K (jMeszar os & Ree^ 

Upon integrating equation , we find that the equi- 
librium tide contribution alone yields a core temperature 



^QlOol/2^1/2 

^ r'cq CQ 



(- 



M 



-1/2 



-5/6 1/2 



V 200a;o 



/. 



100 Hz 



-1/3 



K. 



(50) 



For comparison, linear dissipation of the equilibrium tide 
in a NS-NS binary only heats the NS core to T « 
10^(/gw/100 Hz)5/6 K (|Lai||1994). 

The energy in the dynamical tide mode is rapidly 
thermalized due to nonlinear interactions. Its contri- 
bution to the core temperature is therefore found by 

solving f/(rdyn) = ki''''"^p£^o, which yields Tjyn - 
10^/i(/a/100 Hz)^/6 K. This is comparable to the linear 
result of iLail (|1994[ ) because at frequencies fa < 100 Hz, 
there is enough time for even linear damping to ther- 
malize all the mode energy before the NS merges. We 
therefore conclude that T^q ^ Tdyn. 

6. SUMMARY AND CONCLUSIONS 

To summarize, we considered the nonresonant interac- 
tion between a parent mode and a pair of (p, g) daughter 
modes. We first showed that if the parent mode am- 
plitude I gal > |Kabc|~^, the system is unstable and the 
daughters grow at a rate Fbc — ^c\K,abc<la\- We then eval- 
uated the p-g coupling coefficient Kabc and showed that 
it can be very large when the radial wavelength of the 
p-mode nearly equals that of the y-mode. 

After calculating Kabc for a NS model, we evaluated the 
stability of the tide in coalescing NS binaries to the PGI 
and the resonant parametric instability. We found that 
the equilibrium and dynamical tides are both unstable 
to the PGI at frequencies /gw ^ 20 Hz and excite daugh- 
ter modes that grow 10'^ times faster than the orbital 
inspiral rate ig^. By contrast, resonant parametric cou- 
pling is either stable (in the case of the equilibrium tide) 
or excites daughters at a rate that is much smaller than 
the PGI growth rate (in the case of the dynamical tide). 
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We then considered the saturation of the PGI in or- 
der to determine how it might influence a coalescing NS 
binary. Rather than attempt to solve for the satura- 
tion, which is a difficult problem, we left the saturation 
energy as a free parameter. We showed that if a daugh- 
ter that is driven by the equilibrium tide saturates near 
its wave breaking amplitude, it induces an orbital phase 
error A(/)eq > 1. Since there are A^oq > 100 unstable 
daughter modes, A</)oq » 1 is a possibility. And be- 
cause A(/)eq oc f~^, most of the contribution to the phase 
error comes at large orbital separations (/gw ^ 50 Hz). 
This may have important implications for ground based 
gravitational wave detections because the early inspiral 
contains an important portion of the signal (jCutler et al.l 
11993V We also found that tidal heating due to the equi- 
librium tide PGI can raise the temperature of the NS 
core to ~ 10^° K by /gw ^ 100 Hz. Such early heat- 
ing of the NS might influence the electromagnetic and 
gravitational wave signature of the inspiral and merger. 
Finally, we found that the dynamical tide PGI does not 
significantly alter the linear theory est i mates of A(/)dyn 
or dynamical tide heating found bv iLail ()1994D . 

The saturation of the PGI is a significant source of 
uncertainty and if, instead, the daughters saturate at 
an amplitude well below their wave breaking amplitude, 
the instability will have only a small effect on the or- 
bit and NS. Another important source of uncertainty 
that we briefly discussed is the influence of a superfiuid 
core. Even though the core of a cold NS is expected 
to be superfluid, for simplicity we treat e d it a s a nor- 
mal fluid. Recently, IGusakov fc Kantoil (j2012D showed 
that a superfluid core supports (/-modes with a buoy- 
ancy frequency N that can be much smaller than that 
of a normal fluid. Since the p-g coupling coefficient 
Kabc oc max(a;f,/aJc) oc N~^, a superfluid core may have 
a larger Kate and might therefore be unstable to the PGI 
at larger orbital separations (the local PGI growth rate 
is independent oiujc and therefore N; see eq. [17]). Even 
so, the properties of g-modes in a superfluid core are sen- 
sitive to the core temperature and equation of state and 
it is not clear to what extent Kabc changes in a superfluid. 
While a realistic treatment of the saturation and super- 
fluidity present signiflcant challenges, they are needed in 
order to accurately assess the importance of the PGI in 
coalescing NS binaries. 



The PGI may be important in systems other than 
coalescing NS binaries. Two examples that we briefly 
mention are carbon/oxygen white dwarf binaries and bi- 
naries involving solar-type stars. Using a model of a 
M — O.GMq white dwarf, we find that a pair of {p,g) 
daughter modes can couple to the equilibrium tide with 

a coupling strength ~ 10'^— 10*. This suggests that 
in equal mass carbon/oxygen white dwarf binaries, the 
equilibrium tide might be unstable to the PGI out to or- 
bital periods P w 15 min. In solar type stars, {N/ujqY 
is much larger than in NSs and white dwarfs. As a re- 
sult, a p-mode equals the wavelength of a g-mode (see 
eq. [T]) only for p-mode frequencies well above the solar 
acoustic cutoff frequency (« 60cJo). Since the PGI re- 
quires p-modes that form standing waves, this suggests 
that solar-type stars are stable to this form of p-g cou- 
pling. 

Although we focused on the coupling of (p, g) daughter 
pairs in a stellar core (where Ncg/r is nearly constant), 
there are other potentially interesting forms of p-g cou- 
pling. For example, we showed that a high-order p-mode 
parent can couple strongly to (5,5) daughter pairs. Such 
coupling, which does not require the p-mode wavelength 
to match the g-mode wavelength, might be important in 
systems where p-modes are driven to large amplitudes, 
e.g., by turbulent convection. Another example is the 
coupling of p-modes to g-modes in the region near a 
radiative-convective boundary. Because a g-mode wave- 
length increases as it approaches a convection zone, at 
some point near the radiative-convective boundary the g- 
mode wavelength becomes long enough that it can equal 
the wavelength of a lower-order p-mode. While we find 
that such coupling is fairly weak near the base of the 
convective envelope of solar-type stars, it may be more 
significant in other types of stars, such as those that pos- 
sess convective cores. 

We thank Eliot Quataert for useful discussions dur- 
ing the development of this work and helpful comments 
that improved the manuscript. We also thank Andrew 
Steiner for providing equation of state and composition 
tables. This work was supported by NSF AST-0908873 
and NASA NNX09AF98G. 



APPENDIX 

A. STABILITY ANALYSIS OF p-g MODE COUPLING 

In order to derive the p-g instability criterion (eq. [5]) and growth rate (eq. [5]), assume that the daughters are 
initially at infinitesimal amplitude and the parent is oscillating harmonically at a finite amplitude qa — |<Za|e~*'^*. 
Neglecting linear damping for the moment, the daughter amplitude equations become 

iib + ulqb = uj^K^^q* (Al) 
iic+ujlqc^ujlK^^qt, (A2) 

where Kbc = iiTe^'"* and K = Kabc\qa\- To evaluate the stability of the daughters, let qb.c = Qb,cS^p{db,ct), where 
Qb,c is a constant complex amplitude and Ob.c — s + i<yb,c with s and ab^c complex and real constants, respectively. If 
we set ab -\- (Jc = w, the harmonic time dependence cancels and we find 

[9!+c.!][ie*r+c.'^]^u;lu'jK\'. (A3) 

The real part of this equation gives 



+ [2abac - {ab - a^f + + lo^] + alal - wlal - uj^al + ujIujUi - \K\^) = 



(A4) 
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and the imaginary part gives 



2 Ul'^fJc ~ Lolob 

S = (TbO-c 



These can be combined to give the following more revealing pair of equations 



uji) ± J{ujI - uj^,)^ + Saba, (w^ + a;2) _ Aa^a^iab - + 4w2^2|^|2 



abac {ijjl-ujl) +2{ab~acf{ujl+ijjl) + {ab~act ~ ujluj^\K\^{ab - acf = 



(A5) 

(A6) 
(A7) 



We can use these two equations and the condition ab + a, ~ lo to solve for the three unknowns (s, ctn and a,)- 
For resonant, self-coupled modes we havewft = = io/2 + A and ab = ac = 1^/2, where the detuning A <^ uj. It 

is straightforward to show that equation (|A6[) then yields the usual resonant parametric instability criterion (in the 

absence of linear damping); namely, _Re(s) > if |iir| > 4A/a; (cf. eq. [7]). 
Foi p-g coupling with ujb 3> and \K\ ~ 1, equation (|A7|) implies {ab — ac)^ ~ and aba, — uj'^{uJc/oJb)'^\K\^ <^ 



Equation (|A6I) then reduces to 



(A8) 



and Re{s) > if \K\ > 1; i.e., if \K.abc\ > \Qa\ ^ the daughters are unstable and grow at a rate T^c = Re{s) :^ uJc\K\ — 
\. Furthermore, cr^ ~ uj^ and a'^ ~ ^^(wc/wb)'', and thus during the exponential growth the p-mode oscillates 
at nearly the forcing frequency cj while the g-mode oscillates at a frequency <C oj^uJc- 
When linear damping is included, the equations that determine s, ab, and a, are ab + a, — to and 

[el + ^bOb + u;l] [ie: f + 7c0: + Wc] = ^'b^l\K\^- (A9) 

The real (imaginary) part of this expression is a polynomial equation of order (s^)- Rather than explicitly solve 
these equations, we use the Routh-Hurwitz theorem to determine an approximate instability criterion. The theorem 
states that there exists a solution with Re{s) > if t he o rder-zero term in the polynomial equation is negative. Applied 
to the order-zero term of the real part of equation (jA9p . we find that a sufficient condition for instability is 



\K\' 



1 > 



abac {abac + -fblc) 



(AlO) 



For quasi -adiab atic modes iOb,c > lb,c, and since ab + ac — to <^ tUb for p-g coupling, we expect the right hand side of 
equation (jAlOl) to be less than one (and typically, much less than one). The instability criterion with damping is then 
approximately \K\ > 1, which is the same as that without damping. We have confirmed this result with numerical 
experiments over a range of parameters. 

When the parent is the equilibrium tide, we can still express the stability criterion as \K\ > 1, where now (see 
WAQB) 

-1 

- ^ > ■ M/,,„xrKr^ ( - 1 . (All) 



^bc 



Here W£m = 47r(2^ + 1) ^y£„i(7r/2, 0) and is the Hansen coefficient for the fc-th harmonic of the orbit, i.e., uj — fcfi. 
For a circular orbit and t — 2, \K\ ~ where the tidal amplitude factor e = {M' /M){R/a)^ . 

B. ALTERNATIVE DERIVATION OF THE EQUILIBRIUM TIDE p-g COUPLING COEFFICIENT 

In § 13.11 we stated that for p-g coupling it is straightforward to derive equation (ITlT) starting from the original, 
compact form for Habc- Here we do so for the particular case of an equilibrium tide parent (mode a) coupled to a p-g 
daughter pair (6, c). In original, compact form (see, e.g., eq. [A24] in WAQB), 



^abi 



1 

^0 



X p 



(ri - 1)' 



dlnp 



V • aV bV ■ c 



(Fi - 1) (aj,6^ V • c + 6f,cJ V • a + c^a^ V • b) - Tic^^a'b' c''^.,jk 



(Bl) 



where the subscript semicolon denotes covariant derivative, Fi is the adiabatic index, and $ is the background gravi- 
tational potential. Because we are primarily interested in short wavelength perturbations or the equilibrium tide, this 
expression does not include the terms arising from perturbed gravity (i.e., unlike the full numerical calculation described 
in the main text, here we make the Cowling approximation). In cartesian coordinates with z pointing in the vertical di- 
rection, we have V-a = for the equilibrium tide, V-b ~ dbz/dz for thep-mode, and V-c ~ Cz/H <C dcz/dz <^ dcx/dz 
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for the g-mode (we assume, for simplicity, that the p- and modes do not vary in the y direction). The dominant 
terms in equation (IBip are therefore 



TT I d'x p 



6|; + (ri-i)c>^v.6 



dx dz dz ^ dx dz dz 



Since 5ph/ p — —dbz/dz, the couphng depends on the compression due to the p-mode times the radial shear dcx/dz of 
the g-mode. The expressions for and Cx are essentially the same as those given by br and Ch in equations ([8]) and 
respectively, as long as the daughters' vertical wavelengths are much shorter than their horizontal wavelengths 

1 /2 

and the scale over which background quantities vary. Specifically, bz = At cos (f>b/uJb, where Ai, = [Eoat/ pSxy] , 

ab — cj^ {J '^J^dz^ ^, and 5*0:1, = J dxdy is a constant horizontal area factor that cancels upon integration. Similarly, 

Cx = Ac cos(/)c/wc, where Ac — [EqUc/ pSxy]^^'^ , and ac = (N/z) (/ Ndlnz) ^. For daughters with nearly equal 
vertical wavelengths, kc kh ujh/cs, and 

1 / a daz klAbAc . , . , 

Kabc ~— I a X V ip— sm (pb sm (pc 

ho J dx ujbiOc 

1 /"j3 ^iP Eq .1/2 dazUJb . , . , 

-TT d X — -- — [abac] siiKpbSmcpc 

Eo J pci Axy ox UJc 

~ / dzlabad^^"^ sin (j)b sin (pc- (B3) 

J ox UJc 

For the equilibrium tide daz/dx ~ az/R and we see that this expression for the equilibrium tide p-g coupling coefficient 
(i.e., agrees with equation (1111) modulo the order unity constants of angular integration. 

C. ACOUSTIC CUTOFF FREQUENCY, MODE DAMPING RATE, AND OVERLAP INTEGRAL 

To estimate the NS acoustic cutoff frequency Wac — Cs/2Hp, we use the iRead et al.l ()2009[ ) parametrization of the 
low density EOS of cold matter (valid for the density range p = [10'^, 10^^] g cm~^; see their Appendix C): 

p{p)^K,p^^, e{p)^pc^ + -^, yl^p) = ElEl^ (CI) 

Pi - 1 ^ e+p 



where e is the energy density, Ucq = \/ dpjdt ~ Cs is the equilibrium sound speed, and the parameters Ki and F,; are 
step-functions of p given in Table II of Read et al. For p < 2 x 10^ g cm~'^, have Cg ~ 4.6 x lO^pJj'^^ cm s^^ and 

^ 210 (^) pr°-29c.o, (C2) 
where p4 = p/lO** g cm~^ and 514 = g/lO^"' cm s" -2. For p2 < [(T^)/106 K] \ ideal gas pressure dominates and 



a..^l400 (^j^^j (C3) 

where p is the mean molecular weight. For simplicity, our calculation ignores the finite shear modulus of the NS crust 
and the density discontinuities therein. Since the energy density of the p-modes with frequency near is much larger 
in the core than in the crust (see Figure [5l), we do not expect these properties of the crust to significantly affect the 
value of Wac- 

The linear damping rate, 7a, of a mode in a cold NS is dominated by viscous dis s ipatio n (see, e.g., 
iReisenegger fc Goldreichlll992f ). In order to calculate 7^, we use the method described in iLail ()1994l ) and assume 
that the viscosity is dominated by electron-electron scattering as given by Lai's equation (8.23)13 For high-order 
modes, we find that for the M = 1.4M0 SLy4 NS model, the damping rates are well fit by the formula 
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7a - 7oT,^^An ^ ) , (C4) 



There appears to be an error m Lai s expressions for the shear u-- i ■ iu j ■ * r j iu _i n a ^^ 

J , J . . . 010 J o ir,\ ■ 1 ■ multiplying the damping rates lor modes with m^u. As a result, 

and bulk damping rates (his equations 8.10 and 8.12 : given his t-jo onc i- 4-4-uj ■ 1 c n 

, „ Ti 1 I 1 , I r . r //I , I iM ///I I iM his Bxpressiou 8.25 overestimates the damping rate ot t: = 2,m = 

dehnition ol r^m, there should not be a tactor oi (£+ |m,|)!/(t: — |m|)! ±2 modes by a factor of 24 
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where Tg — T/10'* K is the core temperature. For p-modes (70,0;) = (2 x lO^'* s^^,2) and for g-modes (70, a) = 
(3x10-^8-1,-2). 

In § 14. 2[ we use the linear overlap integral laim = {MR^)~^ f d^xpa ■ V(r^l£m) (see WAQB) to calculate the 

dynamical tide amplitude Igi'^'"^''!. For high-order £ = 2 modes, we find Ia£m — lai'^a/^a)^ , where for p-modes 
(Jo,^) = (10,-3) and for g-modes (/o,/3) = (0.3,2). 

D. LOCAL GROWTH RATE OF THE p-g INSTABILITY 

In order to derive the local growth rate of the PGI, we write the total Lagrangian displacement as a sum of the 
linear and nonlinear displacement ^ = + ^j^j. We can then obtain the equations of motion for (see WAQB): 

pLi = fl [U + 2/2 [Clin, U = fl [U + 2/2 Kun, Cl] , (Dl) 

where fi[^] is the linear force and /2[C,C] is the leading-order nonlinear force, including the nonlinear tide term 
P(Cni ■ V)Vf7. The second equality, which follows because ^ is real, allows us to use WAQB's definition for Kate', 
namely, Kabc = 2-E(7^ / d'^x^^ ■ f2[ibT^c] absorb the factor of 2 multiplying into our definition of the p-g 
coupling coefficient). Since we are interested in the growth rate of the PGI, we neglect the nonlinear inhomogeneous 
term /2Kiin,€iin] (see WAQB for a discussion of this term) and the saturation term /2[^nb^ni] (which is negligible 
when the daughter amplitudes are small). At small daughter amplitudes, the parent determines the linear displacement 
^lin = la^a the (6, c) = {p, d) daughter pair determine the nonlinear displacement ~ qf,^^ -I- (jcCc- Therefore 

+ Uc\=qbfi[ib] + Qcfiiu + 2<z:%*/2 [c^hi + 2<z:g:/2 [cci 

= -p^Ub^b - P^llcic + 2g*g6/2 [C Cb] + 2'7*g*/2 [la, i*c] , (D2) 
where we used the eigenvalue equation fi[£,},] = —pujIS,^,. Defining the radially-local inner product 

{ib,ic)^ I d?^pCb-ic. (D3) 

Jr Jn 

where Ar is many wavelengths long but ^ r, we find from equations ([5]) and © that for short wavelength daughters 
^bi^b^ib) — EoUbAr, i^ciicT^c) — ^ottcAr, and (|;,,Cc) — (the p-mode and g-mode are locally orthogonal even if 
their wavelengths match since, e.g., br oc cost/)?, and Cr oc sin 0c). Taking the local inner product of the equation of 
motion, we obtain the local amplitude equations for the daughters 

dK* 

q, + ulq,^ula-'^q:q: (D4) 
q, + Lolq, = ula-^'^qlql , (D5) 

where we neglected the self-coupling terms because they are small relative to the p-g coupling terms and, since Kabc is 
a slowly varying function of r in the coupling region, we took 



hAr 

73„ 



Ar Jq drdil dr 



From the stability analysis of Appendix we infer that when unstable, the local amplitude equations (jD4[) and (jD5 
yield a local daughter driving rate 

dHabc 



r 1-1/2 
1 6c — i^c [abac] 

The global driving rate is related to the local driving rate by 



<la- 



dr 



(D7) 



/ 1 /2 '"^ 

Ebc ^ l^aKafjcl ^ / dr [a^ac] Ti,c- (D8) 



Since ab^c — Vb \l I ^b l^^' where Vb — Cg and Vc — Lu'^r/KcN are the radial group velocities of the modes, this states 
that the global driving rate is the average of the local driving rate weighted by the geometric mean of the time the 
daughters spend at each radius. 
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